Supersymmetry approach to nuclear-spin-polarization-induced quantum dot structure 
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In nuclear-spin-polarization-induced quantum dots the electrons are confined through local nu- 
clear spin polarization. The model electron confinement potential is time-dependent due to the 
nuclear spin difi'usion and relaxation processes. It can be well-approximated by a Gaussian curve 
which is not an exactly solvable potential. We demonstrate that it can also be approximated by 
multisoliton potentials for the zero value of the angular momentum and by their singular ana- 
logues for other values of momentum without any loss of calculational accuracy. We obtain these 
potentials by supersymmetric (or equivalently Darboux) transformations from the zero potential. 
The main advantage of such potentials is that they are exactly solvable. Time-dependence of the 
nuclear-spin-polarization-induced quantum dot energy levels is found. 

PACS numbers: 73.23.Ab, 72.25.Ab, 75.40.Gb 
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I. INTRODUCTION 

In recent times theoretical and experimental investi- 
gations of quantum dots is attracting a considerable at- 
tention In particu- 
lar, it has been suggested that a novel class of quantum 
dots 0, ^ , so-called nuclear-spin-polarization- induced 
(NSPI) quantum dots (QD), as well as other NSPI low- 
dimensional electron structures 0, 0, , could be cre- 
ated through locally polarized nuclear spins. The main 
idea of creating NSPI structures consists of the follow- 
ing. Electron and nuclear spins interact via the contact 
hyperfine interaction that can be described by an effec- 
tive hyperfine field B/i/, which acts at spins of electrons 
and contributes to the electronic Hamiltonian through a 
Zeeman-type potential g*iiBO'^hf (r, t) /2 Ig^. (Here g* 
is the effective electron gr- factor, /is is the Bohr magneton 
and a is the Pauli matrix- vector). The Zeeman splitting 
results in the potential that is attractive for the electrons 
with one spin projection and repulsive for others. In this 
model the energy is shifted by a constant of the order of 
the Fermi energy by means of a gate. Then the potential, 
created by the inhomogeneous nuclear spin polarization 



Uconf{r,t) 
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liBBhf (r, t) 
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is confining. It is assumed that Bhf (r,t) > 0. 

The idea to create NSPI QDs has been proposed for 
the first time by Fleurov et al |0. They used a pertur- 
bation theory to find the modification of the energy spec- 
trum of a traditional quantum dot due to the polarized 
nuclear spins. Formalism describing NSPI structures was 
thereafter developed and used in investigations of NSPI 



quantum wires . It was established that properties of 
NSPI structures are time-dependent because of nuclear 
spin diffusion and relaxation processes. Subsequently, 
this formalism has been used to investigate NSPI quan- 
tum dots jlljj , NSPI quantum rings 0| , and NSPI peri- 
odic structures Moreover, related low-dimensional 
electron structures created through modulation of spin- 
orbit interaction coupling constant have also been con- 
sidered [13. 

It was assumed in ll'l that a NSPI QD is created in 
the region of the intersection of the two-dimensional elec- 
tron gas with cylindrically polarized nuclear spins. The 
electron energy levels at any particular time were found 
as a solution of the 2D radial Shrodinger equation with a 
Gaussian confining potential. The 2D radial Shrodinger 
equation (as well as 3D one) with a Gaussian potential 
does not give analytic solutions. Different methods were 
used to solve this problem ll7|,|13, 19, 20, 2l| . A parabolic 
approximation of the Gaussian potential has been used in 
Ref. pdl |. However, such an approximation is acceptable 
only for the ground state pdUl^. 

In this paper we extend the study of NSPI systems 
by considering NSPI QD created through spherically po- 
larized nuclear spins. Time dependence of the electron 
confining potential of Gaussian type is found as a so- 
lution of the diffusion equation with a relaxation term 
and the problem is thus reduced to solving the 3D ra- 
dial Shrodinger equation where the time is involved as 
a parameter. For the zero value of the angular momen- 
tum we approximate the Gaussian potential, which is not 
exactly solvable, by a multisoliton potential, which is ex- 
actly solvable. The multisoliton potential is obtained by 
the technique of supersymmetric quantum mechanics (or 
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equivalently by the method of Darboux transformations) . 
For other angular momenta the method gives their singu- 
lar analogues. It should be mentioned that the applica- 
bility of our approach is not restricted only to NSPI QDs. 
It could be used for describing traditional quantum dots 
[iTj l also. 

The paper is organized as follows. We introduce the 
basic equations describing NSPI QD in Sec. In Sec. 
mil we present the method of calculations. Time depen- 
dence of the electron states in NSPI QD is studied in Sec. 
IIVI The conclusions are drawn in Sec. 



II. BASIC EQUATIONS 

Let us consider a semiconductor structure with locally 
polarized nuclear spins. There are two main mechanisms 
leading to the time dependence of the hyperfine field de- 
scribing the nuclear spin polarization: nuclear spin relax- 
ation and nuclear spin diffusion. Then the evolution of 
the hyperfine field is governed by the diffusion equation 



dt 



(2) 



accounting also for the relaxation processes Here 
D is the spin-diffusion coefficient, A is the usual three- 
dimensional Laplacian, and Ti is nuclear spin relaxation 
time 122. In the simplest case, we can assume a 

Gaussian form B^f (r, 0) = Bq exp {^—-^^ for the initial 

condition. The parameters d and Bq define the half- 
width and the amplitude of the initial distribution of the 
hyperfine field respectively. Then the solution of Eq. 
is: 



Bhfir,t)^Boe " l 



(3) 



where to = The maximum nuclear field in GaAs can 
be as high as Bhf = 5. ST in the limit that all nuclear 
spins are fully polarized [S^. This high level of nuclear 
spin polarization has been achieved experimentally. For 
example, the optical pumping of nuclear spins in 2DEG 
has demonstrated the nuclear spin polarization of the 
order of 90% [23|- A similar high polarization has been 
created by quantum Hall edge states (85%) p5]. These 
techniques are developed to allow local polarization and 
controllability of nuclear spins 26| . 

Our approach is based on the following electronic 
Hamiltonian: 

h"^ \o*\ 

where m* is the electron effective mass. The time scale 
introduced by the nuclear spins is several orders of mag- 
nitude larger than the time scale of typical electron equi- 
libration processes. In such a case the electrons feel a 



quasi-constant average nuclear field. This simplifies cal- 
culation by avoiding the complications which would ap- 
pear when solving the Schrodinger equation with the time 
dependence due to polarized nuclei. We take into ac- 
count only the electrons with the spins along the z axis, 
for which the effective potential is attractive. 

The electron energy levels in NSPI QD are determined 
by the radial Schrodinger equation that in appropriate 
units takes the form 



-^ + y(')(x,t)-i. 



-0„,i(a;,t) = 0, 



(5) 



where ipnj {x, t) is related with the radial part of the wave 
function Rn,i(x,t) as follows: ipn.iix^t) = xRn,i{x,t), 



V'^^\x,t) 



l{l + l) Bhf{x,t) 
Bhf (0,0) ' 
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(6) 



X = r/d, 7 = \g*\fiBB,,f{0,0)/ {2Eo), E = e/Eo is the 
energy in dimensionless units, e is the energy, Eq — 
/ (2m* d^), and / = 0, 1, . . .. In the next Section we 
describe the method we are using to approximate the 
potential © with the effective hyperfine field B given by 
Eq. @. 



III. METHOD OF CALCULATION 

The Darboux transformation method also known as 
the method of supersymmetric quantum mechanics is an 
effective tool for solving different problems of theoretical 
and mathematical physics (for reviews see E^). Here we 
are using its property to approximate an interaction be- 
tween composite particles by a local potential with an 
experimental accuracy, which previously proved to be 
very efficient in describing iso-phase (also known as phase 
equivalent) potentials in nuclear physics 2S„^J. Below 
we outline shortly the main features of the method we 
need. 

Suppose one knows all solutions of the Schrodinger 
equation with a given reference potential Vq 



hoii = E^, 



(7) 



Then, using a simple algorithm, one can construct a 
huge multiparameter family of exactly solvable poten- 
tials. The solutions Lp of the Schrodinger equation with 
these potentials are obtained by acting on solutions of 
with a differential operator, ip — Ltjj. In the simplest 
case this is a first order operator 



L 



-d/dx + iu{x) 



(8) 



where the real function w{x) (known in supersymmetric 
quantum mechanics as superpotential) is defined as the 
logarithmic derivative of a solution u of {Tj). So, one has 



w = u' (x) / u{x) , hou — au, 



(9) 
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with a < Eq, where Eq is the ground state energy of ho 
if it has a discrete spectrum or the lower bound of the 
continuous spectrum. The function u is cahed transfor- 
mation or factorization function and a its factorization 
constant or factorization energy. The potential Vi of the 
Hamiltonian hi = —cP/dx^ + Vi, hiip = Eip, is defined 
in terms of the superpotential w as 

Vi{x) ^Vo{x) ~2w'{x) . (10) 

Eq. ^ defines a first order Darboux transformation. In 
the following we shall deal with chains of N successive 
transformations of this type. 

Here we will use very special chains introduced in . 
They are generated by the following system of transfor- 
mation functions 

Vi{x), . . . , Vl{x), Ul + i{x), Vl+l{x), . . . , Un{x), Vn{x) (11) 

hoUj{x) = -a^jU-j{x) , h()Vj{x) = -b'^Vj{x) , (12) 

where Vj are regular {vj(0) = 0) and Uj irregular ('Uj(O) ^ 
0) at the origin. They have arbitrary eigenvalues —aj"^ 
and —bj^ respectively, but always below Eq. If we are in- 
terested in the final action of the chain only, the solution 
ipN{x,k) of the transformed equation with the Hamilto- 
nian 

hN = -(f/dx^ + Vn (13) 

corresponding to the energy E ~ k"^ is given by [30l | 

tpNix,k) = W{ui, . . . ,UN,ipoix, k)) W^^{ui, . . . ,un) 

(14) 

where W are Wronskians expressed in terms of Wj, de- 
noting symbolically any function of (|ll|l and of '0o(2;, k) 
which is a solution of the original Schrodinger equation 
corresponding to the same energy E, N — 2n ~ I. In the 
Hamiltonian H13|l the transformed potential is 

Vn ^Vo-2-^lnWiui,...,UN). (15) 

For = 1 one has W{ui) = ui and one recovers H1U|) 
with u ~ ui. If Vo is finite at the origin, Vn behaves as 
l{l + l)x~'^ when a; ^ 0. Therefore the parameter I can 
be associated with the value of the angular momentum 
in The formulas H14|l and H15|l result from the re- 
placement of a chain of TV first order transformations by 
a single iVth order transformation, which happens to be 
more efficient in practical calculations. 

In Ref. we obtained that the transformed Jost 

function Fn is related to the initial Jost function Fq by 

J = l ^ j = l+l 

For I — the first product is unity. Since a Jost function 
is analytic in the upper half of the complex /c-plane (see 
e. g. [31|) all 6's must be positive whereas the a's can 
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FIG. 1: Comparison between different multisoliton potentials. 
1: two soliton, 2: four soliton, 3: six soliton, 4: eight and ten 
soliton, bold dashed line represent the Gaussian potential. 



have any sign, so that every positive Oj corresponds to a 
discrete level E — —a^ of hN- 

In our case we choose Vq = so that only exponentials 
are involved in the final Hamiltonian. Moreover, since 
the potential ^ is symmetric with respect to inversion 
X — > —x we take Uj{x) — cosh(aja:) and Vj — sinh(6ja:). 
In this case for ^ = we get symmetric multisoliton po- 
tentials (see e.g. the first of Refs. [13 )■ Moreover, the 
Wronskian from (I15|l can be expressed as a sum of hy- 
perbolic cosines 32]. For Z > we get singular at the 
origin analogues of multisoliton potentials. To compare 
them with the second term of the right hand side of © 
we have to subtract the centrifugal part and consider 
Veff — Vn — l{l + l)x^^. The potential Vefj depends 
on N = 2n — I parameters Oj and bj which we fit to the 
Gaussian curve of the second term of the right hand side 
of ijnji for all values of t. This gives us the time depen- 
dence of the parameters aj and bj and hence the energy 
levels Ej = — a^. To get an idea how many functions 
of type Hll|l could be involved in a particular case, we 
plotted the potentials for TV = 2, 4, . . . , 10, Z = 0, 
Bq = l,to/Ti = 1, T = in Fig. 1. It is clearly seen from 
this figure the solitonic structure of the curves for N < 8 
(curves 1, 2 and 3). The result for = 8 (curve 4) is not 
distinguishable of that for = 10 and practically coin- 
cides with the given Gaussian potential curve. Moreover, 
for TV = 10 the highest energy level is equal to zero with 
a high precision. This indicates that this potential has 
only four discrete levels which are not located near the 
ionization threshold E = and possibly one level near 
this threshold. Another remark worth making is that the 
long distance behavior of our potential is ~ exp(— Apcc) 
where Aq = "^j + whereas the potential l|njl tends 
to zero much faster, as ^ exp(— Aix^). 
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FIG. 2: (Color online) Energy spectra of electrons in NSPI 
QD with initial half- width d — 1/im and Bhf{r = 0, t = 0) — 
2.65T as a function of time, Ti/to = 1, Z = (a) and I = 1 
(b). The solid lines are the energy levels obtained using the 
multisoliton approximate potential. The dots correspond to 
the energy levels, that were obtained as a numerical solution 
of the Shroodinger equation with the Gaussian potential. 



IV. RESULTS 

The time dependence of the confining hyperfine field 
given by Eq. Q determines the time-dependence of the 
electron energy levels in the NSPI QD. There are two 
characteristic times in the problem: the diffusion char- 
acteristic time to and the relaxation characteristic time 
Ti. For t ^ to Ti we distinguish diffusive regime, the 
times t such that t ^ to ^ Ti correspond to intermedi- 
ate regime, and when t ^ Ti <^ to we are in relaxation 
regime. Here t is the observation time. 

We found that the time-dependence of the energy levels 
in all regimes are qualitatively similar for NSPI QD with 
spherically symmetric confining potential: the number of 
energy levels, as well as their depth, decreases with time. 
It is interesting to note that the number of energy lev- 
els in cylindrically symmetrical NSPI QD in the diffusive 
regime remains constant and the number of trans- 
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FIG. 3; (Color online) Energy spectra of electrons in NSPI 
QDs with spherical and cylindrical symmetry in the diffusive 
regime, Ti/to = 100. The parameters of calculations are as 
on Fig. H 



verse energy levels in NSPI QW in the diffusive regime 
even increases with time j^. This remarkable fact is 
related to the symmetry of the confining potential: an 
increase of the symmetry results in a faster spreading of 
the confining potential due to the nuclear spin diffusion 
process. 

In order to estimate the accuracy of our approxi- 
mation, we compared the results obtained by super- 
symmetric transformations with numerical solutions of 
the Schrodinger equation with the Gaussian potential. 
We have used the Shooting Method to solve the ra- 
dial Schrodinger equation, subjecting the solution to the 
following boundary conditions: Rn.i (r — s- 0) = r' and 
Rn.i — > oo) — > which correspond to a discrete spec- 
trum eigenfunction. It is convenient to define the func- 
tion Pn,i{r) as Rn,i{r) = 7''P„j(r) thus reducing the first 
of the above conditions to {dPn,i/dr)\r=o = 0. 

Figure 2 shows a representative result of our calcu- 
lations. The time-dependence of the energy levels was 
calculated for NSPI QD based on GaAs with initial half- 
width d = 1/zm and Bhf{r = 0,t = 0) = 2.65T, corre- 
sponding to 50% nuclear spin polarization, in the inter- 
mediate regime, when Ti/io = 1- The number of energy 
levels in the NSPI QD decreases with time. A compari- 
son between analytical and numerical results shows that 
the method of supersymmetric transformations gives a 
good approximation for almost all energies except for a 
tiny interval near E ^ 0. This is clearly seen from Fig. 
2(a) for n = 3 and Fig. 2(b) for n = 2 energy levels. 
Such a difference occurs because of different long distance 
asymptotics of the Gaussian and multisoliton potentials, 
discussed at the end of the Section Ulll 

Faster spreading of nuclear spin polarization in the case 
of 3D diffusion than in the case of 2D diffusion results in 
qualitatively different behaviour of the energy spectrum 



5 



in NSPI QDs with spherical and cyhndrical symmetry in 
the diffusive regime, as demonstrated in Figure 3 for n — 
0. It is seen from Fig. 3 (this was exactly proved in Ref. 

that the number of energy levels in the cylindrically 
symmetrical NSPI QD in the diffusive regime is constant 
in time. In the case of spherically symmetric NSPI QD, 
the number of energy levels decreases in time. 

V. CONCLUSIONS 

We have investigated the electronic structure of a 
quantum dot created through a spherically symmetric 
local nuclear spin polarization within multisoliton poten- 
tial approximation obtained by means of Darboux trans- 
formations. In particular, the electron energy spectrum 
of the NSPI QD having spherical symmetry was calcu- 
lated as a function of time. We found a specific feature 
of the evolution of such NSPI QDs - the number of en- 
ergy levels, as well as their depth, decreases with time 
in all regimes. This observation contrast with previously 
studied cylindrically symmetrical NSPI QD. 

Moreover, it has been demonstrated that the approxi- 
mation of a Gaussian potential by a multisoliton poten- 
tial is much more efficient in comparison with the approx- 
imation by a parabolic potential. The parabolic approx- 



imation is acceptable only for the ground state, whereas 
multisoliton potentials give an excellent agreement with 
the exact results for all energies except for a tiny inter- 
val near E — 0. Moreover, for an electron confined by a 
multisoliton potential the wave function is given analyt- 
ically. We would like also to point out that the method 
presented here can be applied for studying not only NSPI 
QDs but it also can be useful for describing, for instance, 
traditional quantum dots. 

Aknowledgments 

We acknowledge useful discussions with Professor V. 
Privman. This research was supported in part by the Na- 
tional Security Agency and Advanced Research and De- 
velopment Activity under Army Research Ofhce contract 
DAAD-19-02-1-0035, and by the National Science Foun- 
dation, grant DMR-0121146. B. F. S. acknowledges par- 
tial support from President Grant of Russia 1743.2003.2, 
the Spanish Ministerio de Education, Cultura y Deporte 
Grant SAB2000-0240, and the Spanish MCYT and Eu- 
ropean FEDER grant BFM2002-03773. B. F. S. is also 
grateful for the hospitality to the Center for Quantum 
Device Technology of Clarkson University in the Autumn 
of 2003 where this work was started. 



[1] Quantum Dots, edited by L. Jacak, P. Hawrylak, and A. [11 
Wojs (Springer- Verlag, Berlin, 1998). 

[2] F. M. Peeters and J. De Boeck in Handbook of nanos- [12 
tructured materials and nanotechnology, edited by H. S. 
Nalwa, Vol. 3 (Academic Press, New York, 1999), p. 345. [13 

[3] M. A. McCord and D. D. Awschalom, Appl. Phys. Lett. 

57, 2153 (1990). [14 

[4] M. L. Leadbeater, S. J. Allen, Jr., F. DeRosa, J. P. Har- 
bison, T. Sands, R. Ramesh, L. T. Florez, and V. G. [15 
Keramidas, J. Appl. Phys. 69, 4689 (1991); K. M. Kr- 
ishnan, Appl. Phys. Lett. 61, 2365 (1992); M. L. Lead- [16 
beater, C. L. Foden, T. M. Burke, J. H. Burroughes, M. 
P. Grimshaw, D. A. Ritchie, L. L. Wang, and M. Pepper, [17 
J. Phys. Condens. Matter 7, L307 (1995). 

[5] S. J. Bending, K. von Klitzing, and K. Ploog, Phys. Rev. [18 
Lett. 65, 1060 (1990). 

[6] F. M. Peeters and A. Matulis, Phys. Rev. B 48, 15166 [19 
(1993); I. S. Ibrahim, V. A. Schweigert, and F. M. [20 
Peeters, Phys. Rev. B 56, 7508 (1997);J. Reijniers, F. M. [21 
Peeters, and A. Matulis, Phys. Rev. B 59, 2817 (1999). [22 

[7] L. Solimany and B. Kramer, Solid State Commun. 96, 

471 (1995).' [23 

[8] H.-S. Sim, K.-H. Ahn, K. J. Chang, G. Ihm, N. Kim 

and S. J. Lee, Phys. Rev. Lett. 80, 1501 (1998); N. Kim, [24 
G. Ihm, H.-S. Sim and T. W. Kang, Phys. Rev. B 63, 
235317 (2001); H.-S. Sim, G. Ihm, N. Kim, and K. J. [25 
Chang, Phys. Rev. Lett. 87, 146601 (2001). 

[9] S. Nagaraja, J.-P. Leburton, and R. M. Martin, Phys. [26 
Rev. B 60, 8759 (1999). 
[10] V. Fleurov, V. A. Ivanov, F. M. Peeters, and I. D. Vagner, 
Physica E 14, 361 (2002). 



Yu. V. Pershin, Phys. Rev. B, 



,cond- 


mat/0305002 


M. 


Valin-Rodrig 


cond- 


mat/0311077 



A. 



Puente, 



press; 



Serra, 



Yu. V. Pershin, S. N. Shevchenko, I. D. Vagner, and P. 
Wyder, Phys. Rev. B 66, 035303 (2002). 
S.N. Shevchenko, Yu. V. Pershin, and I. D. Vagner, sub- 
mitted to Physica E; cond-mat/0309390 

B. I. Lembrikov, P. Malits, A. S. Kaplunovsky, M. 
Haridim, and I. D. Vagner, to be published. 

C. P. Slichter, Principles of Magnetic Resonance 
(Springer- Verlag, Berlin, 1990), 3d ed. 

J. Adamowski, M. Sobkowicz, B. Szafran, and S. 
Bednarek, Phys. Rev. B 62, 4234 (2000). 
N. Bessis, G. Bessis, and B. Joulakian, J. Phys. A 15, 
3679 (1982). 

C. S. Lai, J. Phys. A 16, L181 (1983). 

R. E. CrandaU, J. Phys. A 16, L395 (1983). 
M. Cohen, J. Phys. A 17, LlOl (1984). 

D. Wolf, Spin- Temperature and Nuclear-Spin Relaxation 
in Matter (Clarendon press, Oxford, 1979). 

D. Paget, G. Lampel, B. Sapoval, and V. I. Safarov, Phys. 
Rev. B 15, 5780 (1977). 

G. Salis, D. D. Awschalom, Y. Ohno and H. Ohno, Phys. 
Rev. B 64, 195304 (2001). 

D. C. Dixon, K. R. Wald, P. L. McEuen, and M. R. 
Melloch, Phys. Rev. B 56, 4743 (1997). 
A. Malinowski, M. A. Brand, R. T. Harley, Physica E 
10, 13 (2001); J.H. Smet, R.A. Deutschmann, F. Ertl, 
W. Wegscheider, G. Abstreiter, and K. von Klitzing, Na- 
ture 415, 281 (2002); T. Machida, T. Yamazaki, and 



6 



S. Komiyama, Appl. Phys. Lett. 80, 4178 (2002); K.R. 
Thurber, L.E. Harrell, R. Fainchtein, and D.D. Smith, 
ibid. 80, 1794 (2002); K. Hashimoto, K. Muraki, T. Saku, 
and Y. Hirayama, Phys. Rev. Lett. 88, 176601 (2002); M. 
Poggio, G.M. Sleeves, R.C. Myers, Y. Kato, A.C. Gos- 
sard, and D.D. Awschalom, cond-mat/0306012 
[27] V. B. Matveev and M. A. Salle, Darboux Transforma- 
tions and Solitons (Berlin: Springer) (1991); F. Cooper, 
A. Khare and U. Sukhatme, Supersymmetry in Quan- 
tum Mechanics (Singapore: World Scientific) (2001); V. 
G. Bagrov and B. F. Samsonov, Phys. Part. Nucl. 28, 
374 (1997). 

[28] J. -M. Sparenberg and D. Baye, Phys. Rev. Lett. 79, 



3802 (1997); J. -M. Sparenberg, Phys. Rev. Lett. 85, 2661 

(2000); H. Leeb, S. A. Sofianos, J. -M. Sparenberg and D. 

Baye, Phys. Rev. C62, 064003 (2000); J. -M. Sparenberg, 

D. Baye and H. Leeb, Phys. Rev. C61, 024605 (2000). 
[29] B. F. Samsonov and Fl. Stancu, Phys. Rev. C 66, 034001 

(2002); B. F. Samsonov and Fl. Stancu, Phys. Rev. C 67, 

054005 (2003). 
[30] M. N. Crum, Quart. J. Math. 6, 121 (1955). 
[31] L. D. Faddeev, Uspehi Mat. Nauk 14, No 4 (88), 57 

(1959). 

[32] B. F. Samsonov, J. Phys. A 33, 591 (2000). 



